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ABSTRACT 



Aims. We present molecular line observations of the massive star forming region IRAS 17233-3606 aimed at studying the molecular 
core associated with the source. 

Methods. The observations were made using the Atacama Pathfinder Experiment telescope in the CO (3-2) and HCO^ (4-3) tran- 
sitions, and in the CH3OH (6^ - 5^), (7^ - 6,^) and CH3CN (16^ - 15^) bands. For the CO(3-2) and HCO+ (4-3) transitions, we 
obtained maps with a size of 70" x 70". The typical angular resolution of the data is ~ 18". 

Results. Our observations reveal an exceptionally rich molecular spectrum, a signpost of hot core activity. Comparisons with two 
other prominent southern hot cores were made through observations in the same frequency setups. We also detected a bipolar outflow 
in CO (3-2) and HCO+ (4-3) lines. Modelling reveals a hot core of size ~ 3" and a temperature of 150 K in the IRAS 17233-3606 
region. The parameters of the molecular outflow are derived through the analysis of the CO (3-2) emission, and are typical of outflows 
driven by high-mass young stellar objects. 

Key words. ISM: jets and outflows ~ ISM: molecules - Stars: individual: IRAS 17233-3606 - Stars: formation 



1. Introduction 

The last decade has seen significant progress in the under- 
standing of how high-mass stars form. Large samples of mas- 
sive young stellar objects (YSOs) were studied with single- 
dish telescopes to investigate their physical properties through 
the analysis of their (sub )mm continuum and molecular emis 



sion ( e.g., Molinari et al]il996ll l998. 2000; Walsh et al.l ll997 , 
19981 [1999,: iHatchell'etain2000u iZhanget al., ,2001.. l2005t 



Sridharanetal.1 l2002t iBeuther et alj l2002bllcl: iFaundez et alj 
2004t IWilliams et all |2004 l2005h . However, an intrinsic fea- 
ture of high-mass stars is that they form in clusters, and that 
most of them are at large (several kpc) distances. Therefore, 
single-dish studies lack the necessary spatial resolution to re- 
solve single protostars. However, they are essential to perform 
statistical studies on large samples of objects, and to identify the 
best candidates for follow up analysis at higher resolution. Since 
the Atacama Large Millimeter Array will regularly operate at 
0.1" scales in the near future, it is now important to select mas- 
sive young stellar object candidates that are observable from the 
southern hemisphere and relatively close-by in order to reach lin- 
ear resolution sufficient for detailed studies of the chemical and 
physical processes in the environment of the newly formed high 
mass stars and the latter's interaction with their environment. 

The region harbouring the prominent far-infrared source 
IRAS 17233-3606 (hereafter IRAS 17233) first came to atten- 
tion through its very intense H2O, OH, and CH3OH masers. 
The combination of the latter is an unambiguous signpost of 
recent hi gh-mass star formation (e.g., iForster & Caswel]||1989l: 
IWalsh et al. 1998) . IRAS 17233, also known as G351.78-0.54, 



Send offprint requests to: S. Leurini 



contains two cm continuum sources with a separation of ~ 12", 
one bei ng a compact HII region, th e other a weak double UCHII 
source (iHughes & MacLeodlll993[) close to the location of the 
maser spots. Recently, higher angular resolution multi-radio 
wavelengths observations by Zapata et al. (2008) found a clus- 
ter of four compact radio objects of sizes equal or less than 0.3" 
and varying (mostly thermal) spectral indices in the region of 
the double radio source. One of them is at the centre of a bipolar 

I II ■ 

outflow traced by OH masers (Fish et al. 2005). Mor eover, non- 
Gauss ian profiles were detected in CO(2-l) bv Osterloh et al.l 
(Il997h which are probably associated with a molecular outflow. 
From these data however, the powering source of the molecular 
outflow remains unidentified. 



Previous studies seem to agree that IRAS 17233 is located 
at the near kinematic distance (betwee n 700 pc and 2.2 kpc, 
iMiettinen et al.ll2006l:lForster & CaswellllT9 89) rather than at the 
far distance (~ 16 kpc, e.g. iMacLeod et aljil99 8). This is sug- 
gested by the high measured intensities of practically all emis- 
sions, which would indicate exceedingly high luminosities if the 
source were at the far distance and, persuasively, by the fact that 
it is at an angular distance of more than 0.5 degrees below the 
Galactic plane. In the following discussion, we adopt a distance 
of 1 kpc fo r the source. Continuu m observations of dust emission 
at 1.2 mm (iFaundez et al.ll2004l) imply a mass of 200 M©, and a 
bolometric luminosity of 2.7 10* L© for the region follows from 
its spectral energy distribution. So far, little data are available 
on thermally excited molecular line emission from IRAS 17233, 
which would allow investigation of the physical and chemical 
status of the region and lay the groundwork for future detailed 
studies with ALMA. 
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In this paper, we report observations of IRAS 17233 with the 
Atacama Pathfinder Experiment telescopq^ in several molecular 
tracers, with typical angular resolutions of 18". These data re- 
veal a bipolar outflow originating fro m a position very c lose to 
the methanol maser spots detected bv lWalsh et al.l (1 19981) . and a 
rich molecular spectrum typical of hot cores near massive young 
stellar objects. 



2. Observations 

The observations were made during 2006 April and August with 
the Atacama Pathfinder Experiment 12 m telescope (APEX) lo- 
cated on Llano de Chajnantor in the Atacama desert of Chile. 
These observations are part of a larger program aimed at study- 
ing ma ssive young stellar objects (YSOs) in the southern hemi- 
sphere (iHieret et al.ll2007[) . 

IR AS 17233 was o bserved with the APEX-2A facility re- 
ceiver JRisacher et al.l |2006) tuned at 356.7 GHz in the up- 
per side band, in order to simultaneously detect the HCO^ 
(4-3) transition (and, from mapping, retrieve the peak posi- 
tion) as well as the CO (3-2) line from the lower side band 
(to detect any outflow activity). A map of 70" x 70" with 
a spacing of 10" was take n centred on t he 1 .2 mm con- 
tinuum peak position from iFaundez et al.l ( |2004|) (a2000 - 
17''26"'43.0', (52000 = -36°09' 15.5"). These observations were 
followed up by ON-OFF integrations on the peak of the HCO^ 
(4-3) emission (aaooo = 17''26'"42.2^ ^2000 = -36°09'21.5") 
in CH3OH (6;^ - 5k, V, = 0,1) at 289.8 GHz in LSB, 
CH3CN (16k - i5K) at 294.45 GHz in USB, and CH3OH 
(7k - 6a:, V, - 0, 1) at 337.8 GHz in LSB. For comparison, 
we also observed the massive star forming regions G327. 3-0.6 
(a20oo = 15''53'«08.6'', (J2000 = -54°37'07.3") and NGC6334(I) 
(a20oo = 17''20"'53.4^ Siom = -35°47'01.5") at 289.8 and 
294.45 GHz. The OFF positions were, in all cases, at (600",0") 
from the reference positions. The APEX facility EFT spectrom- 
eters (Klein etal. 2006) were used in series with 100 MHz of 
overlap, thereby providing 1.8 GHz of bandwidth. Due to the 
high dynamic range of the EFTS (48 dB of the 8 bit analog-to- 
digital converter) the IF processor had to limit the band edges 
with steep band pass filters to avoid aliasing by out-of-band sig- 
nals. However, because of the limited steepness of any analog 
filter, residual aliasing is unavoidable, thereby limiting the effec- 
tive bandwidth of the spectrometer. The region of the spectrum 
which can be affected by this problem is limited to the outer 
-50 MHz of bands. 

In August 2006, continuum cross-sca ns were performed with 
the APEX-2A and the FLASH receivers dHevminck et al.ll2006l) 
at 344.7, 461 and 806.7 GHz respectively, towards the 1.2 mm 
continuum peak position from Faiindez et al. (2004). 

The pointing of the telescope was checked on Jupiter, 
NGC6334(I) and SGRB2-N and found to be accurate to 2"(rms). 

The observed transitions and basic observational parameters 
are summarised in Tables [1] and |2] We used a main-beam efli- 
ciency of 0.73, 0.60 and 0.43 for APEX-2a, FLASH-460 and 
FLASH-810 respectively, to convert antenna temperatures into 
main-beam temperatures. A detailed de scription of APEX and 
its performance is given bv lGtisten et al.l (I2OO61) . 



Table 1. Summary of the parameters of the APEX line observa- 
tions. 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



Line 


Tuning Frequency 


Beam 


Av 


^ sys 




(MHz) 


(") 


(kms-') 


(K) 


HCO+'' (4-3) 


356734.3 


18 


0.05 


270 


CH3OH (6^ - 5^) 


289125.0 


22 


0.06 


293 


CH3CN(16;,-15^) 


294450.0 


21 


0.06 


330 


CH30H(7;,-6;,) 


337800.0 


19 


0.05 


160 



with CO (3-2) in the lower side band. 



Table 2. Summary of the parameters of the APEX continuum 
observations. 



Tuning Frequency Beam 
(MHz) (") 



T 

-* sys 

(K) 



344770.0 
461040.8 
806651.8 



18 
14 



150 

590 

2171 



3. Observational results 

3.1. HCO+ and CO maps 



Figure |l(a)| shows the spectrum of the CO (3-2) transition to- 
ward the peak position of the HCO^ (4-3) emission. The line, 
while at the edge of the observed band, is only marginally af- 
fected by the aliasing out-of-band problem discussed in section 
§ |2] since the 50 MHz at the edge of the EFTS unit starts at a 
velocity of -51 km s ' where the emission in the blue-shifted 
wing is negligible. The line has non-Gaussian blue- and red- 
shifted wings with respect to the systemic ambient cloud veloc- 
ity (-3.4 km s"'), which extend to a maximum blue- and red- 
shifted velocity of ~ +50 km s"' respectively. In addition, two 
narrower lines are detected at ~ 40 and 58 km s ', which could 
be high velocity bullets associated with the outflow, but could 
also be due to unidentified molecular species and be associated 
with the central source. The rest frequencies of these features 
are 345.75 and 345.72 GHz respectively, if they come from the 
lower side band of the receiver, 356.25 and 356.28 GHz if they 
are from the upper side band. 

Broad non-Gaussian wings are also seen in the HCO^ tran- 
sition (vjjj™ = -38, Vj™" = 26 km s"'). However, the profiles of 
the two transitions differ significantly from each other, with the 
CO line showing deep absorption features probably due to self- 
absorption as well as to absorption from other unrelated fore- 
ground clouds. Figures |2(a)| and |2(b)| show a map of the inte- 
grated intensity of the two transitions in the line wings. 

3.2. Molecular observations 

To investigate the physics of the molecular core in IRAS 17233, 
we performed long integration observations toward the peak po- 
sition of the HCO^ emission in several CH3OH transitions and 
in the CH3CN (16a: - 15a:) band, since both molecules can 
be used as pr obes of the physics of the gas (Olmi et al. 19931 
iLeuriniet all 12004). The flexibility of the EFTS at APEX al- 
lows the simultaneous observation of 1.8 GHz bandwidth, with 
a velocity resolution down to ~ 0.05 km s"' at ~1 mm. For the 
analysis of the data, we smoothed all spectra to a resolution of 
1 km s"'. In this way, we observed transitions from the first tor- 
sionally excited levels of methanol together with transitions from 
its ground state; lines from CHj^CN and '''CH3OH also fall in 
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Fig. 1. CO (3-2) (a) and HCO+ (4-3) (b) spectra towards the peak position of tlie HCO+ (4-3) line (marked by a triangle in Fig.|2]i. 
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Fig. 2. a): Map of the integrated intensity of the CO (3-2) transition (blue-shifted emission v - [-50, -10] km s"'; red-shifted 
emission v - [9,50] km s"'; levels from 70 K km s"' in steps of 35, equivalent to ~ 5 cr). b): Map of the integrated intensity 
of the HCO'^ (4-3) transition (blue-shifted emission v = [-35,-8] km s"'; red-shifted emission v = [5,24] km s"'; levels from 
20 K km s"' in steps of 10, equivalent to ~ 5 cr). In both maps dashed contours are used for the blue-shifted emission, solid contours 
for the red-shifted emission. The triangle marks the peak position of the HCO^ emission, the square the compact HI I region, the 
hexagon the centrqid of the multiple radio continuum source, and the asterisk the position of the CH3OH maser spot (IWalsh et alj 
ll998l:IZapata et al]|2008i) . The beam of the APEX telescope at the observing frequency is shown in the left bottom corner. 



the observed bands. We could detect the CH3CN (16a: - 15a:) 
band up to fc - 7. The identification of higher k transitions in 
the observed band is rendered more difficult by their overlap 
with other lines. Having access to high excitation lines, like the 
Vt = 1 CH3OH transitions, and to low opacity lines (e.g., the 
CH^^CN and the CH^^OH bands) guarantees a better determina- 
tion of the physical parameters of the emitting gas (see § 14. 3t . 



The molecular spectrum of IRAS 17233 shows strong emis- 
sion from several CH3OH and CH3CN lines, typical of hot cores; 
moreover, several other transitions from complex molecules are 
detected. For comparison, we observed th e southe rn hot cores 
in NGC6334(I) and G327.3-0.6 (e.g., Sch ilke et al.l l2006) in the 
same setup of IRAS 17233 in the CH3OH (6/r - 5k) and CH3CN 
(16k - 15 A") bands. The spectra are in general comparable, al- 
though the linewidths of IRAS 17233 (Av ~ 6 - 9 km s"') are 
larger than in the other two hot cores (Av ~ 3 - 5 km s"'). 
Several species in IRAS 17233 (CS, H2CO, SO, SO2, HCO+, 
CO) show non-Gaussian wings probably associated with the out- 



flow. Figure [3]presents one of the observed spectral window to- 
wards the three sources. 



4. Physical parameters 

The submm continuum flux of IRAS 17233 is strong enough 
to do continuum cross-scans on the source. These observations 
found the peak emission at (-6", -2") from the position used as 
reference, and not coincident with the peak of the HCO^ (4-3) 
emission (see Fig.|2]i which we used for the spectroscopic obser- 
vations of the hot core. Fit results from averaged cross-scans at 
the diff'erent frequencies are given in Table [3] 

The peak position of the HCO^ emission was derived by 
integrating the line intensity over a range of velocities close 
to the systemic velocity, and contamination from the molecu- 
lar outflow could affect this estimate. The peak of the submm- 
continuum emission is very close (-3". 7, 0". 2) to the p ositio n 
of the methanol maser spots detected by iWalsh et aL J 19981). 
and a weak radio continuum source (named 2b by IZapata et al.l 
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Fig. 3. Spectrum at 294.5 GHz of G327.3-0.6 (bottom), NGC6334(I) (centre) and IRAS 17233 (top). The features labelled in the 
spectra are the ones discussed in section § 14.41 The strong feature at the right end of the bandwidth is a mirror from the lower side 
band of the H2CO (4i,4 - 3 1,3) transition at 281.5 GHz. 



I2OO8I) . With the expected pointing uncertainties the three posi- 
tions are practically indistinguishable. In the following discus- 
sion, we will use the position of the the class II methanol maser 
spots (aaooo = 17''26'"42.6', (^2000 = -36°09'17.7") as ref- 
erence position for the hot core in IRAS 17233, since class II 
methanol masers are found only in the vicinity of massive pro- 
tostars. We also assume that the hot core is the powering source 
of the molecular outflow and that it i s associated with the UCHI I 
region detected at cm wavelengths dHughes & MacLeodlll993h . 
Furthermore, in the following discussion we will assume that 
the total luminosity of the region is driven by the same source 
(harbouring the hot core and the UCHII). The spectral type of 
the UCHII derived from its radio continuum emission is Bl 
dZapata et alj|2008h . while that of the HII region B0.5. However, 
the 850 jum dust continuum map, taken during the ASTROGAL 
survey of the Galactic Plane with the APEX telescope (Schuller, 
priv. comm.), is peaked on the hot core and not on the HII re- 
gion, suggesting that the IR luminosity of the source is indeed 
dominated by the UCHII region. 



Observat ions at shorter wa vele ngths were performe d by sev- 
eral authors: iFix et alj (|1982|) and IWalsh et al.1 (Il999h detected 
emission in the direction of the UCHII region in the near infrared 
K and L bands (2.2jum and 3.3jum respectively), while t he HII re- 
gion is detected in the IR K and L bands and at 10 //m dPix et al] 
119821: IWalsh et alJll999l: lDe Buizer et al."2000^. The area is also 
covered by the GLIMPSE (Benjamin et al. 2003) survey from 
the Spitzer Space Telescope, whose point source catalogue has 
an astrometric accuracy of 0.3". No emission is detected in the 
3.6, 4.5 or 8.0 //m bands towards the hot core, but all bands are 
dominated by strong emission from a source to the north-east of 
the hot core, not associated with any cm continuum emission. 
Emission in the 3.6 and 4.5 fj.m bands is also found close to the 
hot core at (-4.8", 0.2'"). The MIPSGAL survey also covers 
this area dCarev et al.ll2005l) . but both bands (24 and 70yum) are 
saturated. However, by comparing the emission at 24yum (and 
70) with that at 8jum, the peak of the intensity in the MIPS bands 
seems to be shifted towards the hot core position as at mm wave- 
lengths. 



4.1. Continuum 

The mass of the core can be derived from the measured flux 
at submm wavelengths, since in this part of the spectrum the 
continuum emission is dominated by thermal dust emission. 
However, the measurements at 461 and 806.7 GHz should be 
regarded as upper limits to the continuum emission since they 
are likely to be contaminated by emission from the CO(4-3) and 
(7-6) lines respectively. Assuming a distance of 1 kpc for the 
source, and a temperature of 70 K, the flux at 347.7 GHz con- 
verts into a ma ss for the core of 107 Mp fo r a dust emissivity 
of 1.8 cm^g-i ("Ossenkopf & Henning| [T99l. or 130 Mq for a 
standard Hildebrand ( 1983) opacity withjS = 1.5. The value ob- 
tained bv JFaundezet al.1 (120041) was 210 Mq, for T^ = 45 K and 
d=800pc. 

In Fig. m we show the spectral energy distribution of 
IRAS 17233 from 5 GHz to the near IR. Given our resolution, 
several of the IRAC sources in the region could coincide with 
the hot core giving rise to the submm emission or none at all. 
Since we assume that the hot core coincides with the CH3OH 
maser, where no IRAC emission is detected, we actually believe 
the latter to be the case. However, we included the fluxes in the 
IRAC bands in Fig.|4] integrated over the APEX beam. T he cm 
fluxes are from the data presented bv IZapata et al.l (l2008h : they 
correspond to the sum of fluxes of their sources VLA2a-d, re- 
stored with a beam of < 1". An important difference between 
the VLA, APEX and IRAS observations is the spatial resolu- 
tion, which is high in the cm region and coarse at the shorter 
wavelengths. From the source sizes derived from the analysis of 
the submm continuum (see Table |3]l, the emission detected with 
APEX is not confined to a compact source of a few arcsec or 
less, as the VLA sources, but it is rather extended, and probably 
contaminated by emission from components other than the hot 
core, e.g., a surrounding envelope. Therefore, one should regard 
the flux measures from APEX as upper limits to the fluxes of the 
compact hot core associated with the cm free-free emission. 



4.2. Outflow 



Figure |2(a)| shows the map of the molecular outflow originat- 
ing from IRAS 17233 traced by the CO(3-2) line; plotted are the 
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Table 3. Summary of the APEX continuum observations. 



Frequency 
(GHz) 



Tmb 
(K) 



Obs. size 

(") 



Dec. size 
(") 



Flux density 

(Jy) 



347.7 
461.0 
806.7 



1.8 

3.6 

10.7 



22.1 
17.9 
14.6 



12.7 
12.0 
12.4 



204 
1214 



.¥ee-free 

Grey body (70 K) 

Black Ixid; (400 K] 




10* lO" 

I' iGHz] 

Fig. 4. Spectral energy distribution of the hot core in 
IRAS 17233. Observational fluxes from the VLA, APEX, IRAS 
and IRAC are marked as triangles, squares, circles and asterisks 
respectively. For illustration, the flux density curve of an ultra- 
compact HII region, of a 70 K grey body, and of a 400 K black 
body are also given. The solid black line represents the combined 
spectrum of the three components. 



intensities integrated over the velocity ranges covered by veloc- 
ities lower and higher than the systematic velocity, i.e, red- and 
blue-shifted emission, respectively. Although the spectrum at the 
central position (see Fig. |l(a)| l has both red and blue-shifted 
emission at high velocities, implying that the inclination angle 
of the flow with respect to the line of sight is small, the inte- 
grated intensity of the CO and HCO^ lines in the wings shows 
a bipolar distribution: these two properties suggest a medium 
value for the inclination angle of the outflow. The flow has a 
north-south orientation, with the blue-shifted emission towards 
the south, and the red-shifted emission to the north. Strong emis- 
sion at 4.5 //m, which is often interpreted as due to outflows (e.g., 
iNoriega-Crespo et al.l 12004). is detected along the direction of 
the molecular outflow (Fig.|5]l. 

From the CO map it is clear that the lobes are barely resolved 
by our observations, and that higher resolution is needed to de- 
termine parameters such as the collimation factor and the open- 
ing and inclination angles of the outflow. The position angle of 
the flow, computed from the peak positions of the CO integrated 
emission in the two lobes, is close to 0°. 

To compute the physical parameters of the outflowing gas we 
follow ed the standard formalism described by I Cabrit & BertoutI 
( Il990l) . The momentum (p), the energy (E), and the character- 
istic time scale (f) can be derived assuming that the outflow is a 
constant velocity flow with true velocity equal to the maximum 
velocity observed. For the estimate of the mass in the lobes, we 
used a constant value for the integrated intensity of the CO line 
equal to its mean value computed on each lobe starting from 
the Scr level, and assumed that this is indeed representative of 
the flow region. To avoid contamination from the core emis- 
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Fig. 5. Spitzer image of the IRAS 17233-3606 region; the colour 
scheme is the following: 3.6/im (blue), 4.5 yum (green) and 8;um 
(red). The blue and red contours represent the integrated inten- 
sity of the CO(3-2) transition in the blue- and red-shifted wings 
respectively. Contours level are as in Fig. |2(a)| The yellow cross 
shows the position of the CH3OH m aser spot, the square the 
compact HII region dWalsh et alJll998l) . 

Table 4. Parameters of the CO molecular outflow. 



M.o, [Mo] 


0.8 


P[Mq km s"'] 


25.3 


£[10-*^ erg] 


8.6 


l[yr] 


6200 


M,o, [Mo yr-'] 


0.005 


F„[Mokms-'yr-'] 


0.04 


Lm [Le] 


48 



sion, the H2 column density was obtained assuming a relative 
abundance of CO to molecular hydrogen of 10"^ and integrat- 
ing between 9 and 50 km s ' for the red-shifted gas, and -50 to 
-10 km s"' for the blue-shifted gas, which results in lower limits 
to the derived masses. Moreover, part of the red lobe emission is 
not covered by our observations. Since similar studies on molec- 
ular outflows from massive young stellar objects found gas in 
the flow at r elatively high temperatures (Beut heret al. 2002a| ; 
jLeurini et alj |2006), we assume that the excitat ion temperature 
of the outflow gas is 50 K. Some authors (e.g.. lBontemps et al] 
1996) apply a multiplication factor of 3.5 for optical depth ef- 
fects and a mean inclination angle of 57°. 3 in computing the out- 
flow parameters. Since the orientation of the outflow cannot be 
determined from our data, we decided not to correct for the incli- 
nation angle of the outflow. Moreover, observations of CO(3-2) 
and '^CO(3-2) towards other massive flows show that the ratio 
of these t wo lines may change considerably from object to ob- 
ject (e.g.. iLeurini et al.l 120061; iGarav et alj|2007l) thus implying 
differences in the optical depth of CO. The values reported in 
Table|4]should be considered as lower limits to the true values as 
ij no corrections for the optical depth or inclination angle were 
applied; 2) the gas at low velocities was not considered in the 

cal culations. 

iBeuther et al.l (l2002cl) compared the physical properties of 
high-mass outflows to those of low-mass flows by studying the 
correlations between the mass and the luminosity of the core as 
function of the properties of the flow. The flow in IRAS 17233 
has physical parameters characteristic of molecular flows driven 
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by high-mass young stellar objects. In particular, the values de- 
rived for the mass entrainment rate Mtot and for the mechani- 
cal force F are in agreement with a relatively low luminosity 
(L ~ 10"^ Lq) for the powering force of the flow. 

4.3. Molecular core 

The analysis of the molecular spectrum of IRAS 1 7233 was car- 
ried o ut with the XCLASS program (discussed bv lComito et alj 
12005 *). which uses an LTE model to produce synthetic spec- 
tra, and compares the m to the observat ions. The m olecular data 
are f rom the CDMS jMiilleret al .11200 lb and JPL dPickett et alj 
11998.) databases. The parameters defining the synthetic spectrum 
are: source size, rotation temperature, column density, velocity 
off'set and the velocity width. In our case, we did not use any off"- 
set in velocity, but adopted the systemic velocity of the source 
(vlsr - 3.46 km s"') for all molecular species. Several veloc- 
ity components, which are supposed to be non-interacting (i.e. 
the intensities add up linearly), can be used. The program can 
also produce double side band synthetic spectra for comparison 
against the observations. This is needed to model our data, since 
their frequency coverage does not permit sideband deconvolu- 
tion. The galactocentric distance of the source was computed 
by [Ch in et a l.i(ll996l). a nd be found equal to 7.8 kpc. Following 
IWilson & RoodI (1 1 994 . the isotopic ratio '^C/'^'C at this dis- 
tance is equal to 66, very close to the standard value of 60 used 
in XCLASS, and well within the errors of this estimate. 

We started our analysis with the modelling of the CH3CN 
spectrum because several transitions are detected at once, pro- 
viding good constraints on the physics of the gas they trace. The 
corresponding CH^^CN band also falls within the bandwidth of 
our observations; however, at the rest frequencies of these lines, 
blending with other molecular species is severe and it is diflicult 
to assess whether some of the lines are detected or not. However, 
even their non-detection would give us useful upper limits on 
the column density of this species. The best fit was found for a 
source size of 2.7", a temperature of 150 K and a column density 
of ~4xl0'^cm '^ 

IWilner et al.l (1 1995b noticed that for luminosities typical of 
massive young stars, the temperature distribution in spherical, 
centrally illuminated dust clouds can be well approximated by 
the relation 



/ r \0.25 „ 0.4 



(1) 



which relates the temperature of the dust (Tj) at a given dis- 
tance (R) to the luminosity (L) of the inner source. Using this 
equation, and assuming local thermal equilibrium between dust 
and gas, we can independentl y derive the typical angular size 
for the hot core in IRAS 17233. lHughes & MacLeodI (1199 3') esti- 
mated an IR luminosity for the source of 1 .5 x 1 0^ Lq, using a dis- 
tance of 2.2 kpc. This estima te decreases to ~ 2.5 x 10"* Lq when 
adopting a dis tance of 1 kpc (iMacLeod et al.|[l998l) . a value sim- 
ilar to that of iFaundez et al.l (|2004|) . who derived a bolometric 
luminosity of 1.4x10"* L© with a distance of 0.8 kpc. 

Assuming a luminosity of ~ 1.4 x 10"* L©, the diameter of 
the region at 150 K is 2700 AU for IRAS 17233 (see discus- 
sion in § 14.41 ). which corresponds to 2.5" at a distance of 1 kpc, 
in agreement with the value derived through the analysis of the 
molecular spectrum. 

For the best fit values, most of the CH3CN lines are opti- 
cally thick; this implies a degeneracy in the model between the 
source size and the temperature. To derive the 3cr confidence 



levels in the parameters, we fixed the source size to the best 
value derived in our analysis, and performed ax^ analysis for the 
other two parameters in the spectral window between 294.0 and 
294.3 GHz, which contains the CH3CN and CHf CN (16^-15^:) 
lines. We found that the 3 cr confidence level for the kinetic 
temperature of the gas and the column density ranges between 
114 and 290 K, and 2x10'^ and 6x10"' cm^^ respectively. We 
performed the same analysis on the torsionally excited lines of 
methanol. Although excitation by collisions is not an efficient 
way to populate these levels, which are p umped through the in- 
frared radiation field dLeurini et al.ll2007l) . their populations are 
easily thermalised at the temperature of the dust by the dust radi- 
ation field. Since it can be assumed that dust and gas are in ther- 
mal equilibrium (iKriigel & Waknslevll 19841) at the high density 
of the inner region around massive protostars, we believe that the 
LTE approach, defined by one temperature for the dust and for 
the gas, is reasonable also for the v, - 1 lines of methanol. The 
minimum of the x^ distribution is found for different values than 
the ones delivered by the CH3CN analysis (T~ 210 K, s~ 1.7", 
N~ 9 X 10'^ cm"^). However, this model offers a good fit for the 
higher energy levels in the band, while the other lines have high 
optical depths and their profiles do not match the observations. 
Therefore, in Table |5] we list the column density obtained for a 
source size and a temperature equal to that obtained from the 
CH3CN analysis, which still fall in the 3 cr confidence level of 
the;^'^ distribution of CH3OH. In Fig.|6] we present an example 
of the quality of the fit for CH3OH v, - 1 and CH3CN lines. 
The other spectral windows, and the synthetic spectrum corre- 
sponding to the best fit model, are shown in Fig. |8]in the Online 
Material section. 

The parameters used in the model for the different molecular 
species are shown in Table |5] The emission of several molecu- 
lar species is compatible with that of a compact source with the 
temperature and source size derived from acetonitrile. These are 
complex molecules, whose abundances are enhanced by thermal 
evaporation of grain mantles, or whose high excitation lines are 
efficiently excited only by infrared pumping. In the remaining 
cases, the source was assumed to be extended, and the column 
density given in Table |5] corresponds to the beam averaged col- 
umn density. This approach is justified for optically thin emis- 
sion, where the source size and the column density are degener- 
ate parameters of the model. However, it should be noticed that 
this approach was used mostly for molecules where only a lim- 
ited number of transitions were available. 



4.4. Comparison with other hot cores 

Comparison between hot cores associated with sources of differ- 
ent luminosities can help to understand whether their chemical 
composition is a function of the mass and luminosity of the cen- 
tral heating source, whether it simply reflects the evolutionary 
state of the source and its initial conditions. 

The assumed luminosity of IRAS 17233 is ~ 1.4 x 
10"* Lq. For NGC6344(I) and for G327.3-0.6, attempts to de- 
rive their bolometric lumin osity were made by ISandeUl (l2000h 
and IWyrowski et al.l (|2006') respectively, who report values of 
2.6x10^ Lo for NGC6334(I), and 5-15x10"* Lg for G327.3-0.6. 
Therefore, IRAS 17233 is the least powerful of the three hot 
cores we observed. On the other hand, it is likely to be the closest 
among these sources (Dcin.o i ~ 2.9 kpc, Z)ngc6334(I) ~ 1.7 kpc, 
Bergman ll992l;lNeckell 19781) . and it could have the largest beam 
filling factor in our observations. This would explain its excep- 
tionally strong molecular spectrum despite its relatively low lu- 
minosity. 



Leurini et al.: Star formation in the IRAS 17233-3606 region 



3.026 10 



Image Frequenqy (MHz) 
3.025 10 3.024 10 3.023 10 



3.022 10° 



K 

5 




5 



1 - 



O OODO CZ& 



2.828 10 



2.891 10 2.892 10 2.893 10 
Rest Frequency (MHz) 
bnage Frequencv (MHz) 



2.894 10 



2.827 10 



2.826 10 



2.825 10 



2.824 10 




2.941 10 



2.942 10 



2.943 10' 
Rest Frequency (MHz) 



2.944 10 



2.945 10 



Fig. 6. Molecular emission from IRAS 17233-3606. The solid line represents the synthetic spectrum obtained for the best fit solu- 
tion. 



However, using Eq.[Tl we found that the size of the region at 
150 K is equal to 2700 AU for IRAS 17233 (~ 2.5"), 15000 AU 
for NGC6334(I) (~ 8"), and, given its luminosity uncertainty, 
between 5200 and 10 000 AU (~ 1 .7" and 3.4", respectively) for 
G327. 3-0.6. According to these results, the exceptionally strong 
molecular spectrum of IRAS 17233 cannot be attributed to a 
beam filling factor larger than for the other sources. This would 
rather suggest that the strength of the lines from IRAS 17233 are 
due to an intrinsic property of the source. However, the source 
sizes derived through Eq. [1] are affected by the uncertainties in 
the distance, which at least in the case of IRAS 17233 is high, 
and in the luminosity, which for all sources is based on relatively 
low resolution data. 



ISchilke et all ( l2006h and lComito et alj ( l2005h studied the hot 
cores in NGC6334(I), G327. 3-0.6 and Orion-KL with the same 
technique used for IRAS 17233, although at different frequen- 
cies. Therefore, the results of these analyses are affected by the 
same systematic errors of the method, and comparison with our 
results should be straightforward. However, the source size used 
in their studies is not estimated from observations but derived 
from the models. Since this parameter is degenerate with the 



temperature for optically thick lines, and with the temperature 
and the column density for optically thin lines, the column den- 
sities derived for the hot cores are strongly dependent on the 
source size used in the modelling. To cancel out this effect, we 
compared relative abundances in the different sources. Given the 
different spectral windows analysed, only CH3CN is observed in 
all hot cores. Therefore, we used CH3CN as reference molecule 
to compute the abundances of the other species. The species as- 
sociated with the hot cores, and in common between different 
sources, are CH3OH and HNCO (in IRAS 17233, NGC6334(I), 
Orion-KL), C2H5CN (in IRAS 17233, Orion-KL, G327.3-0.6), 
and CH3 OCHO (in IRAS 17 233, NGC6334(I), G327.3-0.6). 
Moreover. lSchoier et alJ (|2002|) studied the molecular content of 
the hot core around the solar type protostar IRAS 16293-2422, 
and derived abundances for four of the species analysed here 
(CH3CN, CH3OH, HNCO and CH3OCHO). 

By comparing the results for these five sources, we could 
not find any trend in the abundances as function of the lumi- 
nosity of the object, and all sources show different values for 
each species. However, for all sources the most abundant species 
is CH3OH, followed in by CH3OCHO, C2H5CN and HNCO. 
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Table 5. Molecular parameters for IRAS 17233. 



Species 


Size 


T,, 


N 


Comment 




(") 


(K) 


(cm-2) 




C2H5CN 


2.7° 


150" 


5(16) 


b 


CH3CCH 


ext." 


100 


2(15) 


c 


CH3CN 


2.7 


150 


4(16) 




CHjOCHO-a 


2.7" 


150" 


2(17) 


b 


CHjOCHO-e 


2.7" 


150" 


2(17) 


b 


CH3OH 


2.7" 


150" 


2(18) 


d 




ext." 


40 


8(15) 




CO 


ext." 


30 


5(19) 


e 


CS 


ext." 


30 


5(15) 


c 


DCN 


ext." 


150 


4(13) 


c 


Ho CO 


ext." 


80 


2(15) 


c 


HoCS 


ext." 


70 


9(14) 


c 


HC3N 


2.7" 


150" 


3(16) 






ext." 


50 


9(15) 




HNCO 


2.7" 


150" 


3(16) 


c 


SO 


ext." 


40 


3(15) 




SO2 


ext." 


60 


5(15) 





Fixed parameter; 

Based on weak or partially blended lines only; 

Based on one or a few lines only; 

Includes vibrationally excited lines; 

Based on C"0 (3-2) only. 



For IRAS 16293-2422, the abundance of CH3OCHO is an up- 
per limit. Figure |7] presents an overview of the molecul ar abun- 
dances with respect to CH3CN in the five sources. Ha tchell et al.l 
(Il998al) performed a similar analysis on a larger sample of hot 
cores, and found no systematic difference between their sources. 
However, we did not include their results in our analysis because 
for most of the molecular species they are limits to the column 
densities. 

We believe that such a comparison of the chemistry of hot 
cores associated with heating objects of different mass and lu- 
minosity, and including hot corinos around low mass protostars, 
would indeed be useful for a better understanding of this evolu- 
tionary stage. However, we stress that this analysis should be 
performed with the same method and on observations of the 
same frequency windows and linear resolutions to cancel out all 
systematic errors. 

Figure [3] presents the molecular spectrum at ~ 294 GHz 
of the three hot cores IRAS 17233, NGC6344(I) and G327.3- 
0.6. The spectra are relatively similar, although NGC6334(I) 
and G327. 3-0.6 show stronger emission than IRAS 17233 in 
the i^CHsOH (6;^ - 5;^) v, = 0, 1 band, while the CH3CN 
{\(iK - 15a-) is more intense in IRAS 17233, at least in its 
lower transitions lines. However, a striking difference between 
the three sources is in the intensity of the SO2 lines (and 
SO transitions at other frequencies), which are by far stronger 
in IRAS 17233 than in the other two cores. The abundance 
of S-bearing molecules can be enhanced by thermal evapo- 
ration of icy grain mantles due to protostellar heating (e.g., 
[Charnley 1997), or by shocks (e.g., Mitchell 1984). Among 
other hot cores, Orion- KL is known to be ex ceptionally strong 
in SO and SO2 (e.g., ISchiIke et allll997h. which a re asso- 
ciated with a molecular outflow JWright et alj Il996h . In the 
case of IRAS 17233, the linewidths of the SO and SO2 transi- 
tions are broad (9-12 km s"'), but reasonably fitted with sin- 
gle Gaussian profiles, suggesting that the bulk of the emission 
does not arise in the outflow. However, for two of the SO2 tran- 
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Fig. 7. Comparison of molecular abundances of CH3OH, 
HNCO, CH3OCHO and C2H5CN respect to CH3CN. 



sitions a two component Gaussian fit gives a better agreement 
with the observations. Both transitions have asymmetric profiles, 
with non-Gaussian blue-shifted emission. Figure [T(b)| shows one 
of the two asymmetric SO2 transitions, which is in the same fre- 
quency setup of the HCO^(4-3) observations, but comes from 
the other side band. The integrated emission of this line is com- 
pact, but slightly elongated along the axis of the molecular out- 
flow. However, only a map of the region in these two molecular 
species at higher angular resolution would unambiguously asso- 
ciate their emission with the hot core or with the outflow. 

THatchefl et al. (1998b) and van der Tak et al. (2003) stud- 
ied several massive star forming regions in different S-bearing 
molecules, and found that the abundances of SO and SO2 
are the ones with the largest spread from source to source. 
Ivan der Tak et all (|2003|) also found that up to 50% of the SO and 
SO2 emission arises in high-velocity gas. No trend in their col- 
umn densities with temperature was found. These studies agree 
on the fact that the abundances of SO and SO2 are strongly de- 
pendent on the source, which suggests that they are related either 
to its evolutionary stage or to its initial conditions as already 
suggested by chemical models (e.g. . Charnlev 1997 for t he de- 
pendence on the evolutionary stage. IWakelam et al]|2004l for the 
dependence on the initial conditions). 



5. Conclusion 

We observed the massive star forming region IRAS 17233-3606 
with the APEX telescope in several molecular transitions of high 
diagnostic value at a wavelength of 1 mm. Our observations re- 
veal a bipolar molecular outflow associated with strong emission 
at 4.5 jim, and an exceptionally rich molecular spectrum from a 
position very close to the centre of the outflow. To better charac- 
terise the nature of this source, we performed continuum cross- 
scans on the source at 0.9, 0.7 and 0.4 mm, and derived the flux 
density of the source over a beam between 18 and 7.6". 

We modelled the molecular emission of the source using a 
LTE approximation, and found that it is compatible with that of 
a hot core of 150 K and a source size of 2.7". The distance of 
the source, and its bolometric luminosity, are not well known. 
However, assuming that the source is at a distance of 1 kpc, and 
that the bulk of the dust is at 70 K, we found that its submil- 
limetre continuum emission corresponds to a mass in the range 
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of 103-240 Mq. We also analysed the molecular emission as- 
sociated with the bipolar outflow, using the CO (3-2) transition 
and assuming that the emission is optically thin and the gas is 
at 50 K. The physical parameters are typical of those driven by 
high-mass young stellar objects. However, follow-up observa- 
tions in high energy and in low optically depth transitions are 
needed for a better characterisation of the outflow. 

If the estimate of the distance to the source is correct, 
IRAS 17233-3606 is probably placed among the less studied hot 
cores around B-type (proto)stars. This seems to be confirmed by 
the parameters of the molecular outflow, which put its powering 
source among the less luminous massive sources. It is therefore 
an interesting target for future high-angular resolution observa- 
tions aimed at investigating the connection between hot cores 
around objects of high mass and hot corinos found around solar 
type protostars. Moreover, given its close distance and its lo- 
cation in sky, IRAS 17233 will easily be accessible in the near 
future for high linear resolution studies with ALMA. 
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3606. The solid line represents the synthetic spectrum obtained for the best fit solu- 
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